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Abstract 

The convergence properties of cycle expanded periodic orbit expressions for the 
spectra of classical and semiclassical time evolution operators have been studied for 
the open three disk billiard. We present evidence that both the classical and the 
semiclassical Selberg zeta function have poles. Applying a Pade approximation on 
the expansions of the full Euler products, as well as on the individual dynamical 
zeta functions in the products, we calculate the leading poles and the zeros of the 
improved expansions with the first few poles removed. The removal of poles tends 
to change the simple linear exponential convergence of the Selberg zeta functions 
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to an exp{— n 3//2 } decay in the classical case and to an exp{— n 2 } decay in the 
semiclassical case. 

The leading poles of the jth dynamical zeta function are found to equal the 
leading zeros of the j + 1th one: However, in contrast to the zeros, which are all 
simple, the poles seem without exception to be double. The poles are therefore in 
general not completely cancelled by zeros, which has earlier been suggested. The 
only complete cancellations occur in the classical Selberg zeta function between the 
poles (double) of the first and the zeros (squared) of the second dynamical zeta 
function. 

Furthermore, we find strong indications that poles are responsible for the pres- 
ence of spurious zeros in periodic orbit quantized spectra and that these spectra 
can be greatly improved by removing the leading poles, e.g. by using the Pade 
technique. 



PACS: 05.45.+b, 03.65.Sq, 02.30.+g 
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1 Introduction 

Trace formulas in chaotic dynamical systems relate phase space averages to sums 
over periodic orbits g, |6|. Exponentiated trace formulas give rise to 

Selberg type zeta functions, named after corresponding expressions arising in studies 
of billiards on surfaces of negative curvature Selberg zeta functions factorize 
further into products of dynamical zeta functions, each one being an infinite product 
over all primitive nonrepeated periodic orbits (ppo) of the system. Finally, the term 
cycle expansion refers to a certain expansion and truncation of (dynamical and 
Selberg) zeta functions into polynomials. Whereas the original trace formulas and 
the infinite products have the same convergence and analyticity properties, cycle 
expanded periodic orbit expressions typically converge much better || f|. 

Calculations can be improved, if the pole structure is known ||. Typically, 
dynamical zeta functions will have poles; for Selberg zeta functions, one can advance 
arguments [||, [H]] that they should be entire and thus ideally suited for numerical 
purposes. We here present quantitative results on the analyticity properties of 
zeta functions for a 2d conservative dynamical system, a point particle elastically 
scattered off three disks placed symmetrically in the plane [11, |l^, [13|, O]. This 



system is ideally suited for such an investigation since (for sufficiently separated 
disks) it is a hyperbolic system with a good symbolic coding (complete binary, once 
the symmetries are factored out). Periodic orbits can conveniently and accurately 
be computed. 

We proceed with a formal definition of the objects investigated: Let p label 
all primitive non repeated periodic orbits, n p their symbolic length, fx p the Maslov 



index (in a billiard fj, p = 2n p ), and J p the linearization perpendicular to the orbit 
with A p the eigenvalue of largest absolute value. For the two degree of freedom 
system considered here J p is a 2 x 2 matrix of determinant 1 so that the other 
eigenvalue is 1/A p . We then consider the classical Selberg zeta function || 



with 

c7 1 w = n( 1 -^ p i A P r 1 v) » ( 2 ) 

p 

and the semiclassical Selberg zeta function [15, [If], 17] 



. 7 m p -irfj.pir/2 | °o 

zw - p t-?S~ i^d-^ } = j5 g " w (3) 

with 

^ rl ( z ) =n( i - znpe_i/i, " r/2 i A pi" l/2 v) - w 

all as functions of z. We avoid use of additional labels distinguishing classical and 
semiclassical zeta functions and hope that it is clear from the context which one is 
meant. Straightforward formal manipulations allow to express each function (l)-(4) 
as a power series C n z n in z, which, when truncated, yields the cycle expansion. 
Where needed, we will abbreviate with t p the contributions of periodic orbits to the 
dynamical zeta function with j = so that C^ 1 = Ilp(l — ^pAp- 7 ). 

The above expressions are correct for maps, the periods of orbits being inte- 
gers. For flows one would replace z Hp by z np exp{iu;Tp} in the classical case or by 
z np exp{iS p (E)/h} in the semiclassical case, expand in a power series in z and con- 
sider the final result as a function of frequency u or energy E, respectively, for z = 1. 
In billiards, the action is given by S p (E)/fr = L p k(E), where L p is the geometrical 



length of the ppo and k(E) = y/2mE/H is the wave number. In addition to Maslov 



phases there can be further phases due to symmetries [18, 19, ^0[. For the case of 
the three disks, the above expressions are correct in the ^-representation; in the 
^-representation there is an additional phase m if n p is odd. 

Periodic orbits for the three disk system have been computed using Newton's 
method on two different maps, one based on direct description (impact parameter, 
scattering cycle) of collisions with the disks and one based on stationarity of action. 
The computations were done for several values of the ratio p = d/R, where d is 
the disk separation and R the disk radius. Symmetry reduced orbits up to sym- 
bolic length 13 have been found in double precision numerics with relative accuracy 
10 — 14 . The exponential of (|l|) and (||) was computed using all orbits and repetitions 
satisfying n p r < N and then expanded in a series C n z n using the recurrence 



relations of Plemelj and Smithies [21, 22 



Fig. 1 shows the results for the classical Selberg zeta function. An apparently 
faster than exponential decay is observed. In contrast, the semiclassical Selberg zeta 
function seems to decay faster than exponential for the first three or four terms, but 
then settles for an exponential decay, see Fig. 2. In the following we will explain 
this difference in behaviour, show its effects in calculations and demonstrate how 
this knowledge can be used to improve calculations. 

We begin in section 2 with a detailed discussion of the convergence of cycle 
expanded zeta functions, including numerical results for the three disk system. In 
section 3 we turn to methods for identification and removal of poles. The effect of 
poles when calculating quantum resonances is discussed in section 4. We conclude 
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with a short summary and comments on other systems in section 5. 



2 Convergence estimates 

Some insight into the behaviour of the cycle expansions of (0) and @ may be 
obtained by considering the special case of a complete binary code with Ao ~ Ai ~ A 
and a factorization of the eigenvalues of the longer periodic orbits A p ~ A^A" 1 , 
where uq and n\ are the numbers of zeros and ones in the symbolic description of 
p. Then the classical dynamical zeta functions take on the form 



^(^-(l-itelAI^A-'j . (5) 
Expanding the product on j in a power series in z, one finds 

Z(z) = U[{r 1 (z)] j+1 ~J2d n z n (6) 



j 



3/2 

with d n ~ A~ n (see below). Similarly, in the semiclassical case one finds 

C7 1 (z)~(l-2z|A|- 1 /2 A ^) (7) 

and 

^)=n<r 1 w~E< aB ( § ) 

with d n ~ A~ n (essentially due to the Euler product formula, see e.g. eq. (89.18.3) 
in Ref. J24|). Because of this rapid decay, these functions are free of poles. 

In the general hyperbolic case, the expansion of dynamical zeta functions does 
not stop with the linear term but rather continues with exponentially decaying 
coefficients, c n ~ f3 n with |/3| < 1. Summing this geometrical series one finds a pole 
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at z = In ref. ||, 10 [3 has been related to the Lyapunov exponents which 

suggests that the poles of G -1 should be compensated by the zeros of Cj+i- To see 
what terms have to compensate in order to provide faster than exponential decay 
of the coefficients in the Selberg type zeta function, let us consider the classical and 
semiclassical cases in some more detail. 



2.1 The semiclassical case 

With each dynamical zeta function expanded in a power series in z, the semiclassical 
Selberg zeta function (||) looks like 

Z(z) = Hcr\z) (9) 

As is often convenient for dynamical zeta functions we distinguish curvature terms 
Ci n of orbits grouped together with shorter shadowing orbits and fundamental contri- 
butions fi of the shortest orbits, which by definition are not approximated by other 
orbits. Superscript labels here indicate the order j of the zeta function. Subscripts 
indicate the size of the terms in powers of lAI" 1 / 2 , where |A| is a typical instability 
of the shortest orbits (obviously, some uniformity in the Lyapunov exponents is 

assumed here). For instance, f[°^ comes from terms of the form | A p | _1 / 2 e _Vp7r//2 

— 1 (i) 
in £ , whence its subscript equals 1. Fundamentals f- with higher j have addi- 
tional powers lApl - -' and thus i = 1 + 2j. The order of the curvature corrections 
cf z n is determined by two factors: the typical size of the terms contributing (about 
|A|- n / 2 -j™) and an additional factor ~ A~ n due to exponential shadowing of long 
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orbits by short approximands, thus i n = (3 + 2j) n. 

To estimate the convergence behaviour, we will now expand (||) in z, Z(z) = 
J2C n z n , and evaluate the leading order contributions in each coefficient C n . One 
sees that there is a considerable difference in order of magnitude between the funda- 
mental term and the first curvature correction, so to begin with one would expect to 
find significant contributions from fundamentals only. In the ideal two scale approx- 
imation A p = A^A™ 1 all curvatures are identically zero, c$ = 0, and the expansion 
looks as follows 

Z(z) = >*) (l-^z) 

= l-(/ 1 (0) + /i 1) + /f + ...), 

+ (fPfP + fPf™ + A (0) /f ) + / 3 (1) / 5 (2) + . • •) z 2 

00 

- (/PW'+zf ..)z 3 + ... = j2 d' n z n • do) 

n=0 

We notice that the leading terms grow in order (I n ) like 1% = 1, I2 = 1 + 3 = 4, 
I3 = 1 + 3 + 5 = 9, . . ., i.e., logd' n ~ I n = S"=o = n 2 , a quadratic exponential 

convergence. 

In the full evaluation of @ we cannot expect the purely fundamental product 
terms to be leading forever. Arranging products like above acording to the sum of 
the lower indices and calling the leading fundamental terms F n 2 we find the following 
expansion coefficients C n : 

C = l 

C 1 = F 1 -fi 1) -fP + 0(7) 

r - p 1 f(°) f( 2 ) M 1 M f (3) , f (i) f (2) jo) nnn , 

8 



r -tp M i fWJo) f(°) fW f(3) i f( 2 ) Jo) _j_ f (o)_(i) , n / 19 ^ 
r> - J°) _u fWJQ) , f (2)Jo) ,(i) ,(2) (o) „ .(o) (i) 

G 4 - -C 12 + /g Cg + /g Cg - /g /g Cg + i<16 + /l C 15 + . . . + C(1SJ 

- J°) j_ fW _(o) , f (2) (o) ,(i),(2) (o) p inc97^ mi 

^5 - ~C 15 + / 3 C 12 + /g C 12 - / 3 J 5 Cg + . . . + ^25 + • • • + <A^7) (.11; 



etc, where 0(n) indicates terms of size ~ A _ra . Up to and including n = 3 the 
convergence is unconditionally quadratic as in the ideal situation. At larger n, 
however, due to the simple exponential decay of the pure curvatures, these as well 
as mixed curvature and fundamental cross product terms have outgrown the pure 
fundamental ones. Unless there now exist efficient additional cancellations within 
complexes of the form 

jo) _ ,(i) fo) (12) 

c 3n J3 c 3n-3 ' \ LZl ) 

raising their order to at least n 2 , a sudden change in the convergence behaviour 
of the semiclassical Selberg zeta function around n = 4 is to be expected. This is 
indeed what is observed in Figure 2. 

2.2 The classical case 

Following the procedure in the semiclassical treatment above we now analyse the 
slightly more complicated classical Selberg zeta function (|l|), 

j'+i 



^)=n[c7 i w] J as) 



Since in the classical case weights are proportional to powers of |A p | , subscripts 
now indicate the size of the terms in powers of |A|~ X rather than |A| -1 / 2 . 
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(?") 

The convergence behaviour in the ideal case c\J = is found as follows: /,From 
a straightforward expansion of 

Z( Z ) = (l-fl%) (l-ff Z ) 2 (l-fP Z f .... £ d nZ - (14) 

n=0 

one obtains the expansion coefficients 
do = l 

d 1 = -ff ) + ... = 0(l) 

d 2 = 2f{ 0) fi 1) + ... = O(3) 

^3 = -/ 1 ( °W ) + ... = 0(5) 

d 4 = 3 / fV 2 (1) / 2 (1) /f + ... = 0(8) 

4 = -3A (0 V 2 (1) / 2 (1) /3 (2) /3 (2) + ... = O(ll) 

4 = / 1 (0) / 2 (1) / 2 (1) /3 (2) / 3 (2 Vf + ... = o(i4) 

d 7 = -4/[ W W W> + • • • = 0(18) (15) 

etc. The growth rule should be obvious: From the jth zeta function (counting Cj" 1 
as the 1st), there are j consecutive contributions to the leading order terms, each 
increasing the order of magnitude by an amount j, i.e., j = j n grows by one over an 
interval of length An = j n and the total growth in order I n is AI n = I n +An—In = jn- 
For large n one thus ends up with the following differential equations: 
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Eq. (|l6|) gives j n = ra 1//2 , which is inserted into Eq. (|i~7|). The solution of the resulting 
equation is the sought for asymptotic relation /„ ~ n 3 / 2 . 
The full evaluation of Eq. (jl^) gives 

C = l 

C X =F X - 2/ 2 (1) - 3/f + 0(4) 

2 = F 3 - 4 0) - / 2 (1) / 2 (1) + 3/f/f + + 4/f/f + 0(6) 

O s = F 5 - + 2/f) C f - 6/{° VP/f + 2/l°) C W + 3/f cf - . . . + 0(8) 

4 = F 8 - 4°) + 2/f 4 0) " / 2 W / 2 (1) cf - 2/f V 2 (1) 4 1] + 3/i 2) 4 0) + • • • + 0(10) 

5 = -4? + 2/«4 0) " / 2 (1) /i 1} 4 0) +F n + ... + 0(12) (18) 

etc. The leading order terms in each C n which have to compensate in order to get 
the faster than exponential convergence n 3 / 2 are now of the form 

42-2/i 1) 41 2 +/i 1) /f ) 41 4 • d9) 

Things work out nicely until order 4; beginning from order 5 additional cancellations 
are needed. As figure 1 shows, these seem to occur in the classical case. 

2.3 Numerical estimate of curvatures 

A rather crude estimate of the individual curvature terms, tpAp , where 

t p = |A p | _1 in the classical and t p = |A p | -1 / 2 in the semiclassical case, was used 
above to obtain the correct order of magnitude for the full curvatures c- , each 
of which being a sum over individuals with the same symbol length. One would 
be able to benefit more from the results of the preceding sections if there were a 
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better estimate of the individual curvatures c p in . Consider again a system 
with binary symbolic dynamics: 

First note that due to uncertainty in the building of complexes like Cqq 011 = 
*ooon - toooih - Moon + Mooi*i and cj$ 101 = tooioi - *ooi*oi it appears necessary 
to collect curvatures with the same number of zeros and ones into a single term 



p 

where the Kronecker <5's select primitive periodic orbits with number of symbols 
no = n and rt\ = m. After this precaution one may make the following ansatz, 



hoping that the essential stability dependence has been correctly extracted, so that 
the prefactors a noni depend only weakly on stability. We aim at finding an approx- 
imation for the prefactors a no m better than a noni ~ 1. 

We have calculated the prefactors a noni up to symbol length 6 for the open three 
disk system, with values of p (= d/R) ranging from 2.5 to 6.0, and with weights 
given by t p = |A P |~ D . The results of the calculations are presented in Table 1 
(D = 1/2, semiclassical case) and Table 2 (D = 1, classical case). The values for 
D = 1 are roughly twice those for D = 1/2 and some variation in parallel with the 
instabilities Ao and Ai are noticable. The dependence on no and n\ seems to be 
roughly binomial with an additional factor no + n\ ; we conclude that the data in 
Tables 1 and 2 should be more or less well approximated by the following formula: 




(20) 




(21) 



riQ+ni-2 



a 



w Dh D (A (p),Ai(p)) (no+m) 



= Dh D (p)B noni 



(22) 



V 



n -l 



J 
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where Kd captures the dependence on the seperation ratio p. Fig. 3 shows the 
rescaled prefactors a noni / DB noni as a function of p, together with the linear fits of 
the data: 

hxiji) « -0.584 + 0.378p (23) 
hi/ 2 (p) « -0.506 + 0.376p . (24) 

The relative deviations from the linear fits are rather large for small values of p, 
but shrink with increasing p; at p = 6 the maximum relative error is less than 10% 
(semiclassical case). The ansatz 7i£>(Ao,Ai) ~ Kd x | Aq Ai | 1//2 with K\ ~ K4/2 ~ 0.2 
is another simple and relatively accurate estimate. 

Eq. ( |22D gives us information with sufficient detail that we may now return to 
the question whether there exist additional cancellations in the classical and semi- 
classical cycle expansions for the open three disk problem. The estimate of the full 
curvatures becomes 



c£> « nDh D E 



^ n-2 ^ 



n =l 



\no \ n—no 

Aq Ai 



\ n - 1 j 



= nDh D X X 1 (X + X 1 ) n - 2 = nDh D X Xif^- 2 , (25) 

where fr> refers to in the classical and in the semiclassical case, respectively. 
For short, we have written Ao = io^o" 1 and Ai = t±Ai . The perhaps surprising 
observation to emerge from Eq. (p5[) is that the leading pole of Cq 1 {z) has to be 
double: 

E c S^ n = Dh D X X 1 f D 1 z E n (fnzf- 1 ~ Z . (26) 

n n I 1 - JZ^j 
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If now (25) is inserted into the semiclassical complexes (|T^), 

" fDcfl, ~ [n-(n-l)] Dh D \ \ifp~ 2 = Dh D \ Q \ x f^ 2 + , (27) 

the terms do not cancel; there still remains a rest of (roughly) the order 0(3n), 
building up a (simple) pole at z p = 1// 3 (1) = (| A 1 _1/2 A^ 1 + 1 A x | - 1 / 2 A^ 1 ) ~ . This 
is in line with our earlier numerical findings that the semiclassical Selberg zeta 
function is not free of poles. However, if we insert the same expression (25) into the 



classical complexes (1£), 



4 0) - + flcfh « [n-2 (n-1) + (n-2)] Dhu^fi-* = , (28) 

the additional cancellations are there. 

Note that the qualitative results above are independent of the choice of weight 
(i.e., the value of D); the different results for the classical and the semiclassical Sel- 
berg zeta function are entirely due to the difference in the power of Cf 1 . The double 
pole of Co" 1 occurs at z p « fp , which is identical to the lowest order approximation 
of the leading zero of ^ (z). By taking the square of ^ l {z) as in the classical 
Selberg zeta function one doubles the leading zero, which then precisely cancels the 
double pole of Co 1 - ^ ne semiclassical case the simple zero of Ci 1 cancels only one 



pole with a simple pole remaining [Eq. (27)]. We study this point further below. 



The conjecture [|| that the position of the poles of the dynamical zeta function 
Q is given by the zeros of Cj+i remains valid, but the order of the poles is not simple 
but double. Returning to the arguments of Artuso et al. ||, [l0| one notes that they 
are rather liberal with the prefactors; and it is precisely in the prefactors that the 
difference between a simple and a double pole resides. 
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3 Identification and removal of poles 



In the case of maps, phase space averages can be related to zeros of zeta functions 
Z( z ) = Y^n=i)C n z n . Practical calculations estimate such zeros from a truncation 
of the series, Fj^(z) = J2n=o C n z n ■ In ideal situations, exponential H Q or even 
faster than exponential convergence || is obtained. The presence of poles destroys 
faster than exponential convergence and makes it more difficult to calculate the 
exact positions of the zeros; consider a simple case where there is one zero and one 
pole: 

F(z) = = (1 - az) (l + bz + b 2 z 2 + ...) 

oo 

= 1 - (a - b) z - (a - b) bz 2 - . . . = C n z n . (29) 

n=0 

We assume that < b < a. A truncation after the linear term gives a value of the 
zero z'q = (a — b)~ . This is obviously a bad estimate of the true value zq = a -1 if 
a and b are of the same order of magnitude. The inclusion of higher order terms 
only slowly improves z' Q , the error being ~ (b/a) N asymptotically. Furthermore, 
additional "ghost" zeros appear: The number of these unwanted zeros equals N — 1 
and they do not vanish to infinity as N grows large - they cluster around the circle 
\z\ = which borders the region of absolute convergence. To see this, consider 
the function 

- nW _ (i-mXi-w") _ (1 _ az) *£ (bz)j (30) 

3=0 

which differs from F^{z) only in the coefficient Cm- It clearly has N — 1 additional 
zeros on the circle \z\ = b^ 1 . 
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For b > a > the situation is still worse; there is no zero of Fn(z) converging 
to zq = a -1 . 

If on the other hand the pole were absent, one would have a polynomial 1 — az, 
which "converges" to its exact form already after the first term in the "expansion" ; 
the zero zo = a -1 is at once correctly determined. By estimating the value of b 
from the asymptotic behavior of the coefficients C n , we could remove the effect of 
the pole: Assume an estimate b close to b with \b — b\ = |e| <C a. Then the function 

(1 - bz)F(z) = 1 - (a - e) z - e (a - b) z 2 \l + bz + b 2 z 2 + . . .1 (31) 

has already in the linear approximation a zero zq = (a — e) 1 close to z$. The 
position of the ghost zeros may be estimated from the function (cf. Eq. (|30|)) 

N 

(1 - bz)F N (z) = (1 - az)(l + (e/b) ; (32) 

3=1 

for large z, it is dominated by the highest power of z, so that its zeros lie on the 
circle \z\ ~ (e/6) _1 / Ar 6 _1 . They tend to infinity as e — > 0. 
Are several poles present, as e.g. in a rational function 

*X*) = S4 , (33) 

Q{z) 

with polynomials P and Q, the zeros \zi\ <C 1^2 1 <C \z$\ <C ... of Q being different 
from the zeros of P, the removal of the leading pole z\ leaves a function where 
Z2 determines the convergence. In the special case P{z) = (1 — az) and Q(z) = 
(1 — b\z) (1 — b%z) with < 62 <C b\ < a, the removal of z\ = b^ 1 pushes the ghost 
zeros to the neighborhood of \z\ = b^ 1 and the error in the linear estimate of zq 
shrinks from b\ [a (a — 6i)] _1 , which may be larger than a -1 , to 62 [o (a — ^2)] 1 *C 
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a -1 . If Z2 in the last example is not much larger than z±, but close in magnitude, the 
improvement is of course only marginal. Thus, if two (or n) leading poles are close 
in magnitude, one has to remove both (all n) before any considerable improvement 
takes place. 

To remove a certain number of poles in a consistent way, one may use a Pade 



approximation [25|: Assume that N + 1 coefficients C n of the function F(z) = 
P{z)/Q{z) = J2n=o^nZ n ar e known and that one wants to know the power series 
expansion of the unknown functions Pm{z) = J2n=o A n z n and Ql(z) = J2n=o B n z n - 
With the normalization Bq = 1, one has to determine M + l + L coefficients A n , B n , 
so that N = M + L is required. The coefficients are determined from Ql( z )Fn(z) = 
Pn-l(z), i.e., 



N L N-L ]0 X < 

£ 0(N-n-n')C n B n/ z n+n = ]T A n z n ; 9{x) = 

n=Q n'=0 n=0 



■ (34) 

1 x > 



By identifying terms of equal power in z, one finds altogether iV + 1 equations for 
theM + L + l = N + l unknowns. 

Once the coefficients are known, the solution of Ql(z p ) = gives approximately 
the position(s) of the dominant pole(s). If the computed value z p converges for 
increasing L (and N), one may feel confident about really having identified a pole. 
Similarly, solutions of Pm{zq) = that converge for increasing L and iV should give 
improved estimates of the zeros of the full function F(z). (Calculating poles is in 
general much more difficult than finding good values for the zeros; even the leading 
pole requires large N.) 

We have applied this to the three disk system. An approximation with linear 
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denominator Q\{z) already gave good results, confirmed and stabilized by compu- 
tations for quadratic and higher order denominators. The positions of the leading 
zeros and poles of the full Selberg zeta function and of the three first dynamical zeta 
functions are listed in Tables 3-6. The data clearly confirm the conclusion drawn 
at the end of the last section that the poles of QJ 1 « coincide with zeros of Cj+i- 
The poles appeared in the computations either as pairs of nearby real values or as 
pairs of complex conjugate values with small imaginary parts; a strong indication 
that they are double. Assuming that they are all real and double, the identification 
of poles could in some cases be done with up to seven digits precision, applying an 
extrapolation scheme on the results given by different orders of approximations N 
and L. 

The improvement in the decay of the coefficients is apparent from Fig. 4, where 
we show the the behaviour of the coefficients A n in the Pn-i(z)/Qi(z) approxima- 
tion of the semiclassical Selberg zeta function (g); in line with the above discussion 
we have approximated the denominator by Qi(z) = 1 — z/zi, where z\ is the lead- 
ing zero of Ci -1 - With the leading pole gone, the faster than exponential decay now 
continues beyond n = 4 out to n = 5 or 6. 

In the classical Selberg zeta function we were not able to identify any pole with 
the Pade technique (N < 13), which would imply that they are all cancelled by zeros 
of higher order dynamical zeta functions. As our numerical results show, however, 
the (leading) poles of the dynamical zeta functions are double; it is evident that 
the double zeros of C-f 2 do cancel the double poles of CJ -1 , as was also shown in 
the preceding section, but the poles of C-f 2 are then quadruple and too many to be 
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completely cancelled by the triple zeros of C^ 3 . Thus we expect the classical Selberg 
zeta function to have poles as well, the leading ones arising not from (q 1 as in the 
semiclassical case, but from Ci ■ Since the magnitude of these poles is rather large, 
we have not been able to extract them directly from a Pade approximation to Z(z). 



4 Real and spurious zeros in quantum spectra 

As mentioned in the introduction most calculations require the zeta functions not as 
functions of z but rather as functions of frequency u and energy E or wavenumber 
k. Consider therefore a case similar to the one above, but with "energy" dependent 
coefficients, C n = C n (E). Put a = 1, b = e~ s , and replace z in Eq. (||) by e iE to 
obtain a function 

l-e iE 

HE) = r—^ J (35) 
it mimics the behaviour of zeta functions in the case of a bounded system. We 
assume 5 > 0, consistent with < b < a above. The function F(E) has zeros 
E n = 27rn, n E N along the real axis and poles E n i = 2Trn' — i5, n € N along the line 
Im(E) = —5, defining the abscissa of absolute convergence in the complex i?-plane. 
We are interested in the properties of a finite order expansion of F and introduce 
therefore a generalized function F(E, z) with the property F(E) = F(E, 1), 

~ . , l-e iE z 

F(E,z) = - 

1 — e w °z 



1 - {l-e~ s )e iE z - (l-e- s )e- s e 2iE z 2 - . . . = ^ C n (E)z n . (36) 

ra=0 



The expansion in z is only formal; at the end one puts z = 1. All results of the 
former section concerning zeros, poles, etc. can now be used, replacing z by e lE , zq 
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by e °, etc., and relating them to the complex S-plane rather than the z-plane. For 
a finite order expansion of F{E) one may therefore immediately state the following: 

(i) Any truncation of the power series puts the zeros of F(E) off the real axis; in 
the linear approximation we obtain e lE = (l — e _<5 ^ , with solutions E' n = E n — ie, 
where e = — log(l — e~ s ); for 5 sufficiently large e « e~ 5 . The imaginary part of the 
energy vanishes with increasing N like e ~ e~ 5N . Large N are required if the poles 
lie close to the real axis. 

(ii) Ghost zeros of F{E) will be found in the neighborhood of the line lm(E) = 
—5: There are N — 1 of them distributed more or less evenly around the circle 
\e iE \ = e 5 ; they satisfy e iE = e^ + ^ with 5j « 6, < fa < 2tt, j = 1, 2, . . . , TV - 1, 
i.e., E = E n j = E n + <j>j — iSj. The ghost zeros converge towards the abscissa of 
absolute convergence and their number goes to infinity with N. 

One may remove poles in the same manner as before; the formal expansion 
F(E, z) =J2 C n (E)z n is approximated with the Pade technique, z is put to 1, and 
left is a power series approximation of the part of F(E) containing the zeros. The 
energy dependent analogue of the two-pole example in the former section demon- 
strates the general tendency: 

(iii) Removing leading poles has the following effect: Ghost zeros are pushed 
down in the negative imaginary direction; main zeros having small imaginary part 
approach the real axis. 

The numerical investigations in this paper were performed for an open hyperbolic 
system, and property (i) has no relevance. Properties (ii)-(iii) are not restricted to 
bounded systems, though. We present numerical evidence that the removal of the 
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leading pole does indeed push what seems to be ghost zeros in the open three disk 
system far down in the negative complex energy plane: Fig. 5a shows the lower 
part of the spectrum of the N = 8 approximation of the semiclassical Selberg 
zeta function together with an improved spectrum where the leading pole has been 
removed (Pade, L = 1). The original spectrum displays, in addition to the main 
zeros which lie relatively close to the real axis, a whole band of lower lying zeros. 
When plotting spectra with higher N, one obtains main zeros at the positions of 
the old ones, but the zeros in the band seem not to stabilize and their number 
increases with N. In the improved spectrum in Fig. 5a a few zeros in the band 
remain with approximately unchanged or even larger imaginary part, while most 
of the others are pushed away; in a certain range around Re(fc) « 13, however, no 
such improvement of the spectrum can be noted. 

To find out whether the remaining zeros correspond to real resonances, we can 
compare Pade improved spectra for different N. Fig. 5b shows improved (L = 1) 
spectra for N = 6 and N = 8. Except for the zeros in the interval Re(k) « 10 — 16 
most remaining zeros in the N = 8 approximation correspond also to a nearby zero 
in the N = 6 approximation and thus seem to be stable. One exception is the zero 
at k « 24.4 — il.5 - we have checked with the iV = 10 approximation though and 
have found a corresponding zero very close to that position. 

To see the correlation between the quantum spectra in Fig. 5a and the position 
of the leading poles, we have plotted the absolute values of the two lowest zeros of 
Ci 1 (E, z) [assuming that they equal the poles of (^(E, z ) an d Z(E, z)] as a function 
of Re(fc) with lm(k) = —2.5 fixed, see Fig. 6. At low values of Re(k) the magnitude 



21 



of the second pole is much larger than the leading one, so that the removal of the 
leading pole has a large effect on the spectrum, which is also observed in Fig. 5a. 
Around Re(/c) ~ 13 the two poles are very close in magnitude, which explains why 
the removal of just one pole only marginally affects the spectrum. For increasing 
Re(/c) the number of leading poles with comparable magnitude increases; it becomes 
more and more difficult to improve the spectrum. 



5 Summary and discussion 

We have presented numerical evidence that the classical and semiclassical Selberg 
zeta functions for the open three disk system have poles. This shows that results 
for 1-d maps |], |2(| cannot be transferred immediately to higher dimensional sys- 
tems, such as the three disk system. In the classical case at least this is somewhat 
surprising, since the system under consideration is an almost ideal hyperbolic 2d 
system (complete binary symbolics, highly unstable periodic orbits with |A p | 3> 1, 
no intermittency) . In the semiclassical case it was clear already from a simple plot of 
the expansion coefficients that the convergence soon settles for a simple exponential 
decay. This was confirmed by 
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(1) a numerical determination of zeros and poles of dynamical zeta functions, 
showing that the leading poles of C~ equal the leading zeros of CJ+i) but 
are double and therefore not completely cancelled; 

(2) a numerical estimate of the curvature corrections of Cq 1 ^ showing that the 
leading pole is double; 

(3) an analysis of the explicit terms in the expansion of Z(z), which with use 
of (2) showed that necessary cancellations between cross terms, in addition 
to those present within the curvatures, do not occur, causing a transition 
in the convergence rate at around n = 4; 

(4) the fact that stable poles could be extracted from Ql(z), the denominator 

of a Pade approximation, Z(z) rj Pn-l(z)/Ql(z), and that the numerator 

Pm-l{ z ) showed improved convergence. 
A plot of the expansion coefficients of the classical Selberg zeta function shows 

seemingly a faster than exponential decay. This is not in contradiction to (2)-(4) 

above as explained at the end of section 2.3: The poles of Q^ 1 are cancelled by 

the zeros of £j~ , but there are not enough zeros from to cancel the poles in 

£j~ . Nevertheless, we have not been able to locate a pole in the classical Selberg 



product by Pade analysis and the tests of Cvitanovic and Rosenqvist |g7| are also 
consistent with an exp{— n 3//2 } scaling. The presence or absence of a pole in the 
classical Selberg zeta function thus remains an open question. 

Since the zeros of dynamical zeta functions are relatively easy to compute, one 
may also quite easily identify the poles, provided there is a 1:1 (here rather 2:1) 
correspondence between poles and zeros of neighbouring dynamical zeta functions. 
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Our numerical findings give evidence that this is most probably the case for the 
leading zeros and poles. With the precision of our investigations, we were not able 
to identify stable next to leading poles from Ql(z) = more than in a few cases and 
then only with 1-2 digits precision; in these cases there did exist a next to leading 
zero within the uncertainty of the pole. 

There are several ways to remove poles. The Pade numerator Pm{z) directly 
gives an approximation of Z(z) or £q (z). One may also first compute leading pole 
estimates z n from zeros of dynamical zeta functions and then multiply J\ n (1 — z/z n ) 
into the expansion of Z(z). Furthermore, the classical Selberg zeta function can be 
used with semiclassical weights to improve the convergence for semiclassical zeros; 
the poles of Cq 1 will then be gone. If, finally, there actually is a m : 1 correspondence 
between the poles of Q~ l and the zeros of Cj+i f° r an h the following construction, 

is free of poles and has zeros which equal the zeros of Co 1 (z) and Z(z). A numerical 
test for the three disk system (m = 2) shows that this zeta function (classical or 
semiclassical weights) really has faster than exponential convergence all the way 
out to the largest n considered (=13). Eq. (|37| ) and related forms require further 
investigations. 

Recently, quantum spectra of bounded chaotic systems have been computed 



from expansions of the semiclassical Selberg zeta function |28| , |29| , |3Cf1 . Despite 
some success, the calculations were made difficult by slow convergence, no clear 
indication that the zeros of the Selberg zeta function approach the real axis as 
./V — > oo, missing quantum levels, and presence of spurious zeros (i.e., zeros of Z(E) 

24 



not associated with exact quantum eigenvalues). The pole induced properties (i)- 
(ii) of the preceding section are here recognizable. We therefore suggest a relation, 
similar to the simple example above, between the position of poles and the distance 
from the real axis of the zeros of Z(E) [or (^^(E)] in bounded systems. In the 
calculations for the anisotropic Kepler problem and the closed three disk system, 
there is one obvious pole, connected with an orbit not realized by the dynamics, but 
for which heteroclinic orbits of arbitrary length exist. This pole has been removed in 
the calculations reported in [29]. But the present investigation suggests that there 
are further poles, not so simply identified. We suspect that they are at least partially 
responsible for the bad convergence of the zeros of Z{E) (in bounded systems there 
are many other sources of trouble, like intermittency and stable islands). Spurious 
zeros present in cycle expanded spectra could be (in many cases at least) nothing 
but the ghost zeros connected to the poles. Improved spectra should therefore be 
obtained by removing the leading poles, as done here for the open three disk system. 
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Figure captions 

Figure 1: Expansion coefficients C n of the classical Selberg zeta func- 
tion (jlj) for different d/R in the open three disk system, (A\- 
representation.) 

Figure 2: Expansion coefficients C n of the semiclassical Selberg zeta func- 
tion (]3|) for different d/R in the open three disk system. (A±- 
representation.) 

Figure 3: Rescaled semiclassical (boxes) and classical (crosses) curvature 
prefactors, a nQni / 'DB noni , for different p (= d/R) in the A\- 
representation. The straight lines are the corresponding lin- 
ear aproximations (Eq. ([24])): hi/ 2 (p) in the semiclassical case 
(broken line), and h\(p) in the classical case (continuous line). 

Figure 4: Improved convergence is achieved in the semiclassical Selberg 
zeta function @ after removal of the leading pole: The fig- 
ure shows the expansion coefficients of (1 — z/z\)Z(z) where 
z\ is the leading zero of Ci ( z )- The data are for the A%- 
representation. 

Figure 5: Semiclassical resonances for the open three disk system for 
d/R = 3, computed using Selberg's zeta function in the A 2 - 
representation. The energy is expressed in terms of k = k(E) = 
\l2mE /%. (a) Spectrum in the N = 8 approximation, without 
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removal (boxes) and with removal (crosses) of the leading pole, 
(b) Spectrum with the leading pole removed, in the N = 8 
(boxes) and N = 6 (crosses) approximation. 
Figure 6: The magnitude of the leading and next to leading zero of the 
semiclassical (i 1 (E,z) as a function of Re(k) with Im(fc) = 
—2.5. The zeros are assumed to equal the leading and next to 
leading pole of Z(E,z). (d/R = 3, ^-representation.) 



30 



Table captions 

Table 1: Semiclassical curvature prefactors a no ,ni of (q 1 (eq. (Q)) for 
the open three disk system (^-representation). 

Table 2: Classical curvature prefactors o: no ^ ni of Co" 1 (eq- ©) for the 
open three disk system (A\ -representation). 

Table 3: Leading zeros and poles of the semiclassical Selberg zeta func- 
tion (|3|) and of the three lowest order dynamical zeta functions 
(||) in the ^-representation, determined from different Pade 
approximations of the respective functions (N < 13). 

Table 4: Zeros and poles of the semiclassical zeta functions in the A2- 
representation. (Cf. Table 3.) 

Table 5: Zeros and poles of the classical zeta functions in the ^-representation. 
Ellipses in column two indicate that no pole could be deter- 
mined from the Pade approximation. (Cf. Table 3.) 

Table 6: Zeros and poles of the classical zeta functions in the ^-representation. 
Ellipses in column two indicate that no pole could be deter- 
mined from the Pade approximation. (Cf. Table 3.) 
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